High-Pressure Phase Diagram in the Manganites: a Two-site Model Study 
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The pressure dependence of the Curie temperature Tc in manganites, recently studied over a wide 
pressure range, is not quantitatively accounted for by the quenching of Jahn- Teller distortions, and 
suggests the occurrence of a new pressure-activated localizing processes. We present a theoretical 
calculation of Tc based on a two-site double-exchange model with electron-phonon coupling inter- 
action and direct superexchange between the t2g core spins. We calculate the pressure dependence 
of Tc and compare it with the experimental phase diagram. Our results describe the experimen- 
tal behavior quite well if a pressure-activated enhancement of the antiferromagnetic superexchange 
interaction is assumed. 
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I. INTRODUCTION 

Rare-earth manganites (Ai-a^A^MnOa, where A is a 
trivalent rare earth and A' a divalent alkali earth) have 
been the object of renewed interest owing to the discov- 
ery of the Colossal Magneto- Resistance (CMR)^ exhib- 
ited by several of these compounds over the 0.2 < x < 0.5 
doping range. The physics underlining the properties 
of CMR-manganites is very rich and not yet completely 
understood)^ although CMR is commonly described in 
the framework of the Double-Exchange (DE) model. 
DE qualitatively accounts for the closely related phenom- 
ena of CMR and temperature-driven transition from the 
paramagnetic insulating phase (high temperature) to the 
ferromagnetic metallic one (low temperature). Neverthe- 
less, the remarkable quantitative disagreement with the 
experimental data indicates that additional and differ- 
ent physical effects must be introduced."* The key role of 
the Jahn- Teller (JT) effect, i.e., a spontaneous distortion 
of the MnOg octahedron, which reduces the energy of 
the system by lifting the degeneracy of the Cg levels, was 
pointed out by Millis et dlA and it is nowadays widely ac- 
cepted. This effect is naturally associated with a sizeable 
electron-phonon (el-ph) coupling, which reduces the elec- 
tron mobility, leading to the formation of small polarons 
and competing with the delocalizing DE mechanism. A 
complete description of manganite properties at micro- 
scopic level is further complicated by the presence of an 
antiferromagnetic (AF) Super-Exchange (SE) interaction 
between the S = 3/2 spins of the localized t2g electrons, 
which is relevant in the electron-doped regime (x > 0.5) 
and could play a role also in the CMR region. 

High pressure is a powerful tool to investigate the 
role of the JT effect in manganites. The application 
of an external pressure produces a symmetrization of 
the MnOe octahedra, thus reducing the lattice distor- 
tions and enhancing carriers mobility. Therefore, ap- 
plied pressure causes an increase of the metal-insulator 
transition temperature Tjm, which in these systems is 



almost coincident with the Curie temperature Tc- An 
almost linear increase of the transition temperature Tc 
with pressure P has been observed in several mangan- 
ites over a moderate pressure range (0-2 GPa)i^ In a 
recent paper the effect of pressure on the metal-insulator 
transition temperature in Lao.75Cao.25Mn03 {x — 0.25) 
has been studied up to 11 GPa by means of optical 
measurements.^ The results clearly show that the linear 
dependence of Tc{P) is limited only to the low pressure 
regime (P < 3 GPa), whereas in the high-pressure regime 
(P > 6 GPa) the Ic(P) curve bends down and it seems 
to approach an asymptotic value. This behavior, consis- 
tently with previous results from high-pressure spectro- 
scopic measurementSj^'^ has been ascribed to the onset 
of a localizing, pressure-activated, mechanism competing 
with the natural pressure-induced charge delocalization. 

The present paper aims at explaining the experimen- 
tal Tc{P) behavior using a simple theoretical approach 
which enables a deep understanding of the microscopic 
mechanisms driving the observed pressure behavior. We 
propose a model which contains, besides the standard DE 
and el-ph coupling, also AF-SE coupling between the t2g 
core spins, which may provide the new localizing mecha- 
nism at high-P invoked to account for the experimental 
behavior. ^'^'^ It has been proposed that the competition 
between AF-SE and DE accounts for some of the mag- 
netic properties of the manganites.— The possible role 
of the AF-SE term in providing the high-pressure extra 
localizing channel is, indeed, quite reasonable, since the 
superexchange between t2g orbitals may be thought as 
proportional to the square of some hopping integral tt 
between the t2g orbitals of neighboring manganese ions. 
The hopping integral, arising from the overlap of the 
atomic wavefunctions is expected to increase with in- 
creasing pressure. 

An extremely simple, yet efficient, model to compute 
Tc{P) is the two Mn-site cluster, which represents the 
minimal model including DE, el-ph coupling, and SE in- 
teractions. The simplicity of this model enables to carry 
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out exact calculations at finite temperature, which are 
fundamental to describe the experimental P — T phase 
diagram. The two-site model (TSM) here presented has 
already been applied in Ref. where it has been shown 
to successfully catch the relevant physics of manganites 
at T = owing to the extremely short-range nature of 
the interactions occurring in these systems. Here we ex- 
tend the study to finite temperature and discuss the de- 
pendence on pressure of the various parameters in the 
Hamiltonian. In this way we treat, at the same level, the 
quenching of the JT interactions and the role of AF-SE 
coupling in determining the pressure dependence of the 
metal-insulator transition. We show that the experimen- 
tal results are well described if the AF-SE interaction is 
assumed to be renormalized proportionally to the square 
of the hopping matrix element. 

The paper is organized as follows: in Sec. ^] we in- 
troduce the model and the details of the calculation; In 
Sec. mil the pressure dependence of Tc is reported and 
discussed, while Sec. II VI is dedicated to the conclusions. 



II. MODEL AND CALCULATIONS 

The TSM has been extensively employed in the study 
of polaronic systems, since it is the simplest model able to 
describe the crossover from a metal to an almost localized 
small polaron, which is reflected in both ground stateiSii^ 
and the spectral properties ii^ii^ The DE-TSM including 
el-ph coupling has been already studied at T = 0^^ and 
its relevance to the case of CMR manganites has been 
discussediiiiiSiii In particular, the correct ground state 
has been obtained in several coupling regimes. The quali- 
tative difference between quantum (5'=3/2) and classical 
(S = oo) spin cases was also pointed out, showing the 
importance of quantum fluctuations of the core-spins in 
a proper study of manganites^ In the present paper, we 
extend the model of Ref. ^ to finite temperature and 
focus our attention on the temperature-driven crossover 
from paramagnetic to ferromagnetic ordering occurring 
in manganites. In such a minimal cluster no phase tran- 
sition with long-range order can occur but, owing to the 
very short-range nature of the interactions in mangan- 
ites, predictions on possible instabilities taking place in 
thermodynamic limit can be inferred. 

The hamiltonian of our TSM reads: 

<y 1=1,2 

+JiSi ■ §2 — g (rii — 7^2) (a + a)) + ujao^a. (1) 

In the electronic part (first three terms), t is the hop- 
ping integral between the Cg levels, c\ ^ {ci,a) is the cre- 
ation (annihilation) operator for an electron of spin a on 
site i, Jh is the Hund's rule coupling, cfi = cl^aapCip 
is the spin operator on site i {(Jap are the Pauli matri- 
ces). Si is the local spin ( Si = 3/2) due to localized t2g 



core electrons on site i, and Ji is the AF-SE coupling. 
The phonon contribution (last two terms) contains an 
Holstein coupling and an harmonic term, where g is the 
electron-phonon coupling, rii = So-cf o-Ci.CT is the electron 
number operator at site i, and (a) creates (annihi- 
lates) an Einstein phonon of frequency wq, coupled to 
the density difference between the two sites. The above 
form of the el-ph coupling is equivalent to local phonons 
coupled to the on-site electron density n^, after eliminat- 
ing the symmetric phonon mode which couples to the 
total density For the sake of simplicity, we replace the 
JT coupling with the standard Holstein el-ph coupling. 
It has been shown that, as far as the evaluation of Tc is 
concerned, the Holstein coupling in a single orbital model 
gives results extremely close to a two-orbital model with 
JT interactions^^ The independent parameters of this 
model are Jn/t, Ji/t, wq/^, and the dimensionless el-ph 
couphng A = 2g'^ /uot. As already mentioned, the ex- 
treme simplicity of the model does not require further 
approximations and allows for an exact solution for ar- 
bitrary values of the parameters. As an example, the 
quantum nature of phonons and core spin can be fully 
taken into account. 

In the present calculation we will consider the case of 
a single electron on the two sites, for which the Hamilto- 
nian has a very small (64 x 64) electronic Hilbert- 
space (including the quantum S = 3/2 spins of the 
t2g electrons), while the infinite phonon Hilbcrt space 
necessarily requires a truncation up to some maximum 
phonon number Uph- Except from the extreme adiabatic 
case {idQ/t <C 1) and/or extreme multiphononic regime 
{Xt/ujQ 3> 1), where a really large number of phonons 
can be excited, convergence of the relevant states can 
be, in general, achieved for relatively small (<100) Uph- 
In the non-adiabatic and intermediate coupling regime a 
complete diagonalization of |^ can be carried out numer- 
ically. In the case of manganites, where we can estimate 
u!o/t «0.5 and A « 1 (see below), the convergence for the 
low energy states is achieved already at 71^/1= 10. 

The complete diagonalization of H, i.e., the determina- 
tion of all eigenvalues En and eigenvectors \tpn), enables 
the calculation of the thermal average of the nearest- 
neighbor spin correlation operator ^ — (Si ■ S2)/ S'^ as 

^(V„lA^IV'«)exp(-i;„/r) 

KT) = r^TTn ■ (2) 

X;exp(~i;„/r) 

n 

This quantity measures the short-range magnetic corre- 
lations and it is positive for ferromagnetic phases, neg- 
ative for antiferromagnetic, and vanishes in the high- 
temperature paramagnetic phase. As we will show in the 
next section, a simple analysis of the temperature depen- 
dence of the spin-correlation ^(T), allows us to give a re- 
liable estimate of Tc . The choice of this "short-range" es- 
timator of Tc is, in our opinion, preferable to thermody- 
namic estimates in our TSM. The reliability of the TSM 
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in the case of manganites can be assessed also through the 
analogy with a completely different theoretical approach, 
the Dynamical Mean Field Theory (DMFT)^ which has 
been applied to a lattice model for the manganites by Mil- 
lis and coworkers.^*' The DMFT is a powerful nonpertur- 
bative approach which freezes spatial fluctuations, but 
fully retains the local (single-site) quantum dynamics, 
and becomes exact in the limit of infinite coordination^ 
The DMFT maps the original lattice model in the ther- 
modynamic limit onto a self-consistent impurity model 
which interacts with a quantum bath. As already dis- 
cussed in Ref. ITll we emphasize the conceptual analogy 
of these two approaches. In the TS model, the quantum 
nature of the problem is also completely retained, and 
each of the two sites "feels" the presence of the other 
site similarly to the way the impurity site feels the bath 
within DMFT. We stress that, within the DMFT ap- 
proach, the system is in the thermodynamic limit. Both 
the methods, despite their differences, are thus expected 
to well describe the physics of short-range correlations 
in manganites. The choice of the TSM allows us to 
easily include the nearest-neighbors SE-AF interaction 
which is instead unaccessible by the single-site DMFT 
and requires an Extended-DMFT study, where the dy- 
namical spin correlation function is also self-consistently 
evaluatedfi^ or a cluster-DMFT approach, which retains 
short-range dynamical correlations, but substantially in- 
creases the computational effortiSi 



III. PRESSURE DEPENDENCE OF Tc 

In order to carry out the comparison between the ex- 
perimental results reported in Ref. and the Tc{P) val- 
ues which can be determined using the TSM, we need 
to know the pressure dependence of the input model 
parameters t{P), Jh{P), ^oiP), A(P), and Ji(P) for 
Lao.75Cao.25Mn03 at least over the 0-11 GPa pressure 
range. In Ref. 22 the hopping integral t and the Hund's 
coupling Jh for Lao.eTCao.ssMnOs (x — 0.33) have been 
derived from first principles calculations as a function 
of a/flQ, where a is the cubic lattice parameter and 
flo its ambient pressure value. It is reasonable to as- 
sume the same dependence of t upon a in the x = 0.33 
and X — 0.25 materials up to an overall factor, namely 
ix=o.25(a/ao) = Ctx=o.33ia/ao). The proportionality fac- 
tor C can be evaluated exploiting the proportionality be- 
tween t and Tc at P = [Ref. i-e. through the rela- 
tion i,=o.25/ix=o.33 = Tc{x = 0.25)/Tc{x = 0.33). The 
pressure dependence of the Hund's coupling Jh can be 
safely assumed to be identical for x = 0.25 and x — 0.33 
compounds, since Jh is much larger than the other en- 
ergy scales, and the results of the calculation are therefore 
only very weakly dependent on its variations. Finally, 
using the lattice parameter a(P) recently measured over 
the 0-15 GPa pressure range for Lao.75Cao.25Mn03r^ the 
a/ao dependence can be simply converted into a pressure 
dependence. The increase of the hopping integral t{P) as 
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FIG. 1: Pressure dependence of the hopping integral (dashed 
line) and of the dimensionless electron-phonon coupling (solid 
line, see text). 



a function of pressure is plotted in Fig. ^in the relevant 
pressure range (about 12 % at 12 GPa). 

The pressure dependence of the JT phonon frequency 
LUo (P) , directly obtained from Raman measurements car- 
ried out over the 0-14 GPa ranged is also useful to get 
an estimate for A(P). For a JT mode, it is reasonable to 
assume that the electron-phonon coupling g is given by 
g = ^/2Mujodt{a) /da, where M is the ionic mass. Us- 
ing the above set of pressure dependent parameters, the 
estimate for the pressure dependence of the el-ph cou- 
pling A — 2g'^ /luqI can be obtained. The resulting A(P) 
is plotted in Fig ^ The zero-pressure value thus ob- 
tained, A(0) = 1.04, can be checked against an indepen- 
dent estimate of g. The typical Jahn- Teller distortion 
xq = (x) = l/y/2MiJo{{a + a^)) can be estimated, in 
the polaronic regime, as xo — \/2/MuJag/LUo. Using the 
value Xo = 0.09 A for Lao.75Cao.25Mn03 at P = [Ref. 
01, one has A(0) = 0.94, which compares well with the 
previous estimate. The remarkable pressure-induced re- 
duction of A(P) (about 50% at 12 GPa), is consistent 
with the observed enhancement of the metallic charac- 
ter of the system and with the increase of the transi- 
tion temperature The last ingredient is the AF-SE 
coupling Ji with its pressure dependence. To our knowl- 
edge, no direct estimate of this parameter, either exper- 
imental or theoretical, is available. An estimate of the 
zero-pressure value of Ji{P) can be obtained by simply 
imposing the calculated Tc{0) be equal to the known ex- 
perimental value of 220 K for Lao.75Cao.25Mn03 at zero 
pressurei^ 

In Fig. El we show n{T) calculated in the TSM us- 
ing the known zero pressure values of the parameters 
(t(0) « 0.20 eV, wo(0)A(0) = 0.39, Ji/(0)/t(0) = 12.8, 
A(0) = 1.04) and a value of Ji(0) = 0M7t which is con- 
sistent with the experimental value of Tc- The crossover 
from a ferromagnetic state at low temperature, with a 
large positive value of /i, to a paramagnetic state at high 
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FIG. 2: Temperature dependence of the spin-correlation 
(thick solid line), calculated for A =1.04, Jh /t =12.8, 
ujo/t =0.39, and Ji/t =0.047. Light solid and dashed lines 
represent the half-maximum and linear extrapolation meth- 
ods for determining Tc, respectively (see text). 



FIG. 3: Tc(P) calculated with the two-site model for Ji = 
const, (solid line), J\ oct (dashed line), and J\ oc (dotted 
line), compared with experimental values (squares). Empty 
circles represent the DMFT estimate— (see text). Dashed- 
dotted line is an exponential fit of experimental datai. 



temperature with /i ^ is evident. The absence of an 
abrupt variation of fi{T) is expected since the real phase 
transition occurring in the thermodynamic limit cannot 
take place in a finite system and it is replaced by a smooth 
crossover. In the TSM the crossover is quite broad and 
the estimate of the critical temperature is not completely 
straightforward and implies some degree of arbitrariness. 
For this reason, we propose, and compare, two differ- 
ent methods for evaluating Tc- The first one, the half- 
maximum method, simply defines the Tc as the tem- 
perature for which /i(Tc) = [^i{0) - At(oo)]/2 = At(0)/2, 
i.e., for which the spin-correlation becomes one half of its 
zero-temperature value. The second method is based on 
a linear extrapolation around the inflection point of the 
/i(T) curve. These two methods are schematically repre- 
sented in Fig. 121 from which it is apparent that the latter 
method provides a higher value for Tc- Interestingly, the 
two estimates of Tc are actually proportional over a wide 
range of parameter values. Therefore, the choice of one 
method basically affects only the absolute value of Tc 
and not the internal comparison between results obtained 
with a given estimator. In the following, we use the half- 
maximum estimate to determine the effect of pressure 
on Tc- It is important to note that the precise Ji(0) 
value is influenced by the method chosen for determining 
Tc- For example, the determination of Tc with the half- 
maximum method leads to Ji(0) = 0.047^(0) whereas 
the linear extrapolation provides Ji(0) « 0.060^(0). The 
slight difference between the zero pressure values is not 
important since we are mostly interested in the pressure 
dependence of the parameters rather than in their abso- 
lute values. 

As mentioned before, since the pressure dependence of 
Ji is unknown, we have carried out calculations oiTc{P) 
for three different pressure-dependences of Ji{P) namely: 
i) Ji(P) = const.] ii) Ji(P) cx i(P); iii) Ji(P) oc t^{P). 



Bearing in mind that t{P) is an increasing function of 
the pressure (see Fig. the three choices correspond to 
reduce, keep constant, and increase the ratio Ji/t upon 
increasing pressure. The results of the TSM calculations 
together with the experimental data from Ref. are 
shown in Fig. |31 Before discussing the accuracy of our 
results with respect to the experiment, we compare with 
previous DMFT results of Ref. |23, in which a DE term 
with classical core spins, a JT-type el-ph interaction with 
classical phonons and no AF-SE term have been used. 
By combining the Tc{X)/t values reported for n = 0.75 
(i.e., X — 0.25)^" with the present estimate for X{P) and 
t{P), the Tc{P) are readily obtained. The Tc{P) val- 
ues, also shown in Fig. O have been divided by a factor 
1.1 in order to force the agreement with experimental 
data at zero pressure. Such a small difference for the 
Tc{0) value confirms the reliability of the present model 
parameter estimates. Comparing the theoretical results 
shown in Fig. |21 it is interesting to note that DMFT and 
the TSM with Ji cx t{P) provide extremely close results. 
Beside the previously discussed relationship between the 
two approaches, this agreement can be explained as fol- 
lows: The pressure dependence in the DMFT approach is 
basically determined by A, while in our TSM A, loq, and 
Ji depend on P. Since LUo/t is only weakly dependent on 
pressure, the choice Ji (x t leads to a TSM in which the 
pressure dependence is substantially determined only by 
A as in the DMFT estimate. This argument explains the 
observed equivalence between the two approaches and, 
above all, it points out how the pressure dependence of 
the el-ph interaction alone does not provide a satisfactory 
description of the observed experimental high-pressure 
behavior. 

The comparison between our theoretical estimates and 
experimental data is reported in Fig. |31 The pressure 
dependence of the AF interaction remarkably affects the 



5 




P (GPa) 



FIG. 4: Empty circles: pressure dependence of Ji extracted 
from the experimental data of Ref. by means of the two- 
site model (see text). Solid line is a guide to eyes. For sake 
of comparison t^{P) (dashed line) and t{P) (dotted line) are 
also shown. All curves are normalized to their zero-pressure 
value. 



calculated P — T phase diagram. For small pressures, 
the cases i) [ Ji = const] and ii) [Ji oc t] nicely follow the 
experimental Tc, while case iii) [Ji cx t^] gives a lower 
critical temperature. By increasing the pressure, in the 
two first cases Tc increases too rapidly with pressure, 
while the third estimate closely follows the experimental 
data up to 11 GPa. In particular, Ji = const produces an 
almost linear behavior of Tc{P) over the whole pressure 
range. The ability of the Ji cx t^ data to reproduce the 
actual experimental behavior is really encouraging, since 
the t^ dependence may be seen as the most "physical" 
approximation for Ji. The basic hypothesis is that the 
pressure dependence of the overlap between the t2g or- 
bitals is similar to the dependence of the same overlap for 
the Cg orbitals, i.e., that for each pressure tt{P) oc t{P). 
Under this assumption, the superexchange between core 
spins Ji (X t1/U (where U is an effective repulsion en- 
ergy depending on Jh, the Mn-0 charge transfer energy, 
and the on-site Coulomb repulsion) becomes in turn pro- 
portional to t^ /U , t being the hopping integral between 
the Cg orbitals. These findings clearly show that the ef- 
fects of pressure are not only limited to a reduction of 
the el-ph coupling but a remarkable pressure dependence 
of the AF interaction has to be introduced to account for 
the experimental behavior. The same information can be 
obtained by extracting an "experimental" dependence of 
the antiferromagnetic term on pressure Ji (P) . This sim- 
ply amounts to determine, for each pressure value, the 
value of Ji which provides the experimental Tc, given 
that all the other parameters are fixed by the above men- 
tioned conditions. 

In Fig. 0] we show Ji{P), normalized to the zero- 
pressure value together with t'^ (P) / 1'^ (0) . The two quan- 
tities compare quite well, strengthening the Ji{P) oc 
t^{P) ansatz. For sake of comparison the quantity 



t{P)/t{0) corresponding to the case ii) is also shown in 
Fig. 0] The above results strengthen the hypothesis that 
the saturation of Tc{P) observed above 6-7 GPa^ can 
be ascribed to the onset of a regime in which the AF 
interaction is no more negligible and competes with the 
DE. This hypothesis is supported by several recent high- 
pressure experiments. The occurrence of AF interactions 
in a DE system leads to charge localization, which in turn 
may favor the coherence of the JT distortions. This pic- 
ture is consistent with high-pressure X-ray diffraction ex- 
periments, in which the onset of a coherent JT distortion 
has been observed at about 7 GPa in Lao.75Cao.25Mn03 
and other compounds of the Lai_a;Caa;Mn03 series i^^i^^ 
It is also worth to notice that the onset of the AF phase 
in LaMnOa is accompanied by a remarkable softening of 
the B2g-JT phonon.^^'^^ Even if in our case no real AF 
ordering occurs, the enhancement of the AF interaction 
suggested by our analysis is consistent with the observed 
saturation of the same phonon in Lao ysCao 25Mn03 at 
about 7 GPaS. 

The Ji(P) oc <^(P) dependence suggests that, upon in- 
creasing the pressure, the system evolves from a regime 
dominated by DE, where the bond compression leads to 
an increase of Tc, to an intermediate regime, where Tc 
is almost constant and independent on pressure, owing 
to the competition between SE and DE, and eventually 
to a very high pressure regime dominated by the SE 
contribution where Tc starts decreasing. This picture 
is absolutely consistent with the pressure dependence of 
Tc{P) recently observed in several manganites where the 
three regimes are apparent as well as with the 

data of Ref. for Lao.75Cao.25Mn03, where the system 
only reaches the intermediate regime for the maximum 
pressure. The extent of the pressure range over which 
Tc is almost pressure i ndep endent is different for the 
samples studied in Refs. I28l 29 30, although the onset of 
saturation takes place around 4 GPa for all these sam- 
ples. This pressure value is lower than that observed in 
Lao.75Cao.25Mn03i and the difference can be ascribed to 
the presence of strong coupling and/or high cation disor- 
der which have been shown to have strong effects on the 
pressure dependence of Tc^ 



IV. CONCLUSION 

We studied the evolution of the transition tempera- 
ture of Lao.75Cao.25Mn03 with pressure by exactly solv- 
ing a two Mn-site model which contains all the relevant 
microscopic interactions acting in the manganites. The 
theoretical results have been compared with the phase 
diagram recently measured over the 0-11 GPa pressure 
range. ^ As input parameters for the model we used pres- 
sure dependent data available in literature, except for 
the direct AF-SE between the t2g spins, for which, in 
absence of experimental estimates, we compared three 
different pressure dependences. We have shown that the 
pressure dependence law for Ji is indeed crucial to ob- 
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tain a good agreement with the experiment. Neglect- 
ing the effect of pressure on the AF-SE term leads in 
fact to a sizeable overestimate of the Curie temperature 
for P > 6 GPa. On the other hand, the theoretical re- 
sults follow the experimental data over the whole pressure 
range for Ji{P) oc t^{P), while the calculation largely 
overestimates the critical temperature at large pressures 
if a weaker pressure dependence of Ji is assumed. 

The present results give a clear indication that the ob- 
served high pressure behavior of Tq can be ascribed to 
a pressure-induced enhancement of the AF-SE which en- 



ters in competition with the DE above 6-7 GPa. The idea 
of a pressure activated AF interaction is, at least, com- 
patible with the previous experimental evidences of the 
onset of cooperative JT effec t ^-^^-^ and phonon softening^ 
in the very high pressure regime. 
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